Canine Circovirus in Foxes from Northern Italy: Where Did It All Begin?

Canine circovirus (CanineCV) is a recently identified virus affecting both domestic and wild carnivores, including foxes, sometimes in presence of severe clinical signs. Its circulation in wild animals can thus represent a potential threat for endangered species conservation and an infection source for dogs. Nevertheless, no data were available on its circulation in the Alps region of Northern Italy. In the present study, samples collected from 186 foxes in the period 2009-2020 from Valle d'Aosta and Veneto regions were tested using a real-time PCR assay, demonstrating a viral circulation of approximatively 2-5%, depending on the considered regions. Two complete or almost complete genome sequences were obtained, highlighting that the detected strains were part of a so defined "fox only" clade, which suggests that, despite common contact opportunities, Alps foxes are not involved in frequent transmission events to domestic dogs. Such genetic isolation could be at least partially attributed to some sort of independent evolution occurred in the foxes, leading to species barrier. Additionally, CanineCV strains in foxes from Italy were unexpectedly related to those previously identified in foxes from the United Kingdom and Scandinavian area. Combining the history of fox distribution in Europe since the last glacial maximum (LGM) with the viral history allowed us to speculate a long-standing coexistence between European canine circovirus and this host, justifying the peculiar geographic distribution and evolutionary paths of the fox infecting clade.

Tracking the molecular evolution and transmission patterns of SARS-CoV-2 lineage B.1.466.2 in Indonesia based on genomic surveillance data.

BACKGROUND: As a new epi-center of COVID-19 in Asia and a densely populated developing country, Indonesia is facing unprecedented challenges in public health. SARS-CoV-2 lineage B.1.466.2 was reported to be an indigenous dominant strain in Indonesia (once second only to the Delta variant). However, it remains unclear how this variant evolved and spread within such an archipelagic nation. METHODS: For statistical description, the spatiotemporal distributions of the B.1.466.2 variant were plotted using the publicly accessible metadata in GISAID. A total of 1302 complete genome sequences of Indonesian B.1.466.2 strains with high coverage were downloaded from the GISAID's EpiCoV database on 28 August 2021. To determine the molecular evolutionary characteristics, we performed a time-scaled phylogenetic analysis using the maximum likelihood algorithm and called the single nucleotide variants taking the Wuhan-Hu-1 sequence as reference. To investigate the spatiotemporal transmission patterns, we estimated two dynamic parameters (effective population size and effective reproduction number) and reconstructed the phylogeography among different islands. RESULTS: As of the end of August 2021, nearly 85% of the global SARS-CoV-2 lineage B.1.466.2 sequences (including the first one) were obtained from Indonesia. This variant was estimated to account for over 50% of Indonesia's daily infections during the period of March-May 2021. The time-scaled phylogeny suggested that SARS-CoV-2 lineage B.1.466.2 circulating in Indonesia might have originated from Java Island in mid-June 2020 and had evolved into two disproportional and distinct sub-lineages. High-frequency non-synonymous mutations were mostly found in the spike and NSP3; the S-D614G/N439K/P681R co-mutations were identified in its larger sub-lineage. The demographic history was inferred to have experienced four phases, with an exponential growth from October 2020 to February 2021. The effective reproduction number was estimated to have reached its peak (11.18) in late December 2020 and dropped to be less than one after early May 2021. The relevant phylogeography showed that Java and Sumatra might successively act as epi-centers and form a stable transmission loop. Additionally, several long-distance transmission links across seas were revealed. CONCLUSIONS: SARS-CoV-2 variants circulating in the tropical archipelago may follow unique patterns of evolution and transmission. Continuous, extensive and targeted genomic surveillance is essential.

Role of Q675H Mutation in Improving SARS-CoV-2 Spike Interaction with the Furin Binding Pocket.

Genotype screening was implemented in Italy and showed a significant prevalence of new SARS-CoV-2 mutants carrying Q675H mutation, near the furin cleavage site of spike protein. Currently, this mutation, which is expressed on different SARS-CoV-2 lineages circulating worldwide, has not been thoughtfully investigated. Therefore, we performed phylogenetic and biocomputational analysis to better understand SARS-CoV-2 Q675H mutants' evolutionary relationships with other circulating lineages and Q675H function in its molecular context. Our studies reveal that Q675H spike mutation is the result of parallel evolution because it arose independently in separate evolutionary clades. In silico data show that the Q675H mutation gives rise to a hydrogen-bonds network in the spike polar region. This results in an optimized directionality of arginine residues involved in interaction of spike with the furin binding pocket, thus improving proteolytic exposure of the viral protein. Furin was predicted to have a greater affinity for Q675H than Q675 substrate conformations. As a consequence, Q675H mutation could confer a fitness advantage to SARS-CoV-2 by promoting a more efficient viral entry. Interestingly, here we have shown that Q675H spike mutation is documented in all the VOCs. This finding highlights that VOCs are still evolving to enhance viral fitness and to adapt to the human host. At the same time, it may suggest Q675H spike mutation involvement in SARS-CoV-2 evolution.

Virological Characterization of the First 2 COVID-19 Patients Diagnosed in Italy: Phylogenetic Analysis, Virus Shedding Profile From Different Body Sites, and Antibody Response Kinetics.

BACKGROUND: The pathogenesis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection remains unclear. We report the detection of viral RNA from different anatomical districts and the antibody profile in the first 2 COVID-19 cases diagnosed in Italy. METHODS: We tested for SARS-CoV-2 RNA clinical samples, either respiratory and nonrespiratory (ie, saliva, serum, urine, vomit, rectal, ocular, cutaneous, and cervico-vaginal swabs), longitudinally collected from both patients throughout the hospitalization. Serological analysis was carried out on serial serum samples to evaluate IgM, IgA, IgG, and neutralizing antibody levels. RESULTS: SARS-CoV-2 RNA was detected since the early phase of illness, lasting over 2 weeks in both upper and lower respiratory tract samples. Virus isolate was obtained from acute respiratory samples, while no infectious virus was rescued from late respiratory samples with low viral RNA load, collected when serum antibodies had been developed. Several other specimens came back positive, including saliva, vomit, rectal, cutaneous, cervico-vaginal, and ocular swabs. IgM, IgA, and IgG were detected within the first week of diagnosis, with IgG appearing earlier and at higher titers. Neutralizing antibodies developed during the second week, reaching high titers 32 days after diagnosis. CONCLUSIONS: Our longitudinal analysis showed that SARS-CoV-2 RNA can be detected in different body samples, which may be associated with broad tropism and different spectra of clinical manifestations and modes of transmission. Profiling antibody response and neutralizing activity can assist in laboratory diagnosis and surveillance actions.